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Therefore, taking tangent function of both sides of (1), we get 

sin p tan a -{- sin p tan (co — a) , 

1 — sin'p tan a tan (co — a) ' 

or, by dropping functional characteristics of a, co and X, 

sin»(a + : i )-*-(l — sin 3 p— ) = A (2) 

*\ 1 + acoJ \ - 1 + acoj v ' 

By expansion and reduction, (2) becomes 

(co sin p — X sin 2 p)a 2 — Xco cos s pa = X — m smp. ... (3) 

The two roots of a in (3) are 

col cos 2 p . [4(X — co sin p)(co sin p — X sin 2 p) -\- AWcos 2 j)]^ -i 

2(<o smp — Asin 2 p) 2(cosinp — Xsitfp) ' 

coX cos 2 p [4(X—co sin p)(co sin p — X sin 2 p)+AWcos 2 p] Ji 

2(«sinp — Xsin 2 p) 2(cosinp — Xsin 2 p) ' 

which shows, as in a plane, that the point D may be posited at the same 
distance from B that it is from A. 

This problem, in its most general form, has a somewhat remarkable his- 
tory. A general solution may be found in Newton's Universal Arithmetic, 
prob. 24 ; and in Gergonne, Annales de Mathematiques, torn. 10, p. 205. 



ON AD J USTMENT FORMULAS. 



BY E. L. DE FOREST. 

When a series of observed numbers shows irregularities resulting from 
errors of observation, the simplest mode of correcting them is by graphical 
construction, plotting the terms of the series on paper as ordinates to a curve, 
drawing a smooth curve through the points thus obtained so as to coincide 
with them as nearly as may seem best, and then measuring the ordinates 
thus corrected. This will often give results accurate enough for practical 
purposes. It also has the advantage of bringing the conditions of the 
problem together before the eye, so that graphical construction will always 
be a valuable adjunct and preliminary to the use of other methods which 
may be expected to give a higher degree of accuracy. 

Of the analytical methods for adjusting a series, the simplest and most 
easily applied by persons of moderate mathematical knowledge is that which 
presupposes that the terms of the given series are equidistant, or if not, 
that they have been reduced to equidistant ones either graphically or by 
simple interpolation, and then adjusts the middle term u of any group of 
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2m-j-l terms by means of the formula 

u' Q = l u +l 1 {u 1 +u_ 1 )+l 2 {u 2 +u_ 2 )+ . . . +l„(u m +u_ m ), . . (1) 
where l , l lt l 2 &c, are numerical coefficients determined in advance. It 
is most advantageous to assume that the true law of the series is algebraic 
and of a degree not higher than the third, so that if there were no errors of 
observation, the fourth and higher orders of differences would be zero. In 
point of fact, the true analytical law is usually unknown, but the above 
assumption respecting it will be practically safe, provided that the form of 
the given series, when plotted on paper is seen to approximate to that of an 
algebraic curve of the third degree, or a lower degree, throughout every 
group of 2m +1 consecutive terms. Such a curve can have only one point 
of inflexion, so that the method is obviously unsuited to the adjustment of 
series which are known to be of sinuous form, so as to have very frequent 
changes of flexure. 

Subject to the condition named, there may be an infinite number of 
systems of values assigned to the coefficients I, such that if the adjustment 
formula is applied to a series of the third or any lower order, the adjusted 
value u'q will be the same as the given value w . Schiaparelli has investi- 
gated some properties of these formulas, with a view to the reduction of 
meteorological observations chiefly, in his treatise Sul modo di ricavare la 
vera espressione delle leggi della natura dalle curve empiriche; Milan, 1867. 
He showed that the coefficients I must be connected by these two relations ; 

1^+2^+32*3+ ... +mH m =0f w 

He also ascertained what the numerical values of I must be, for values of 
2m+l from 5 to 11 inclusive, in order that the probable errors of the ad- 
justed terms may be reduced to a minimum. The results thus reached are 
really identical with those obtained by the method of least squares, when 
the 2m +1 terms are each supposed to satisfy approximately the equation 

u = a + 6* + ex 2 + cfer 1 
and are also regarded as of equal weight, that is, equally liable to error. 
(Compare an article by the present writer, in the Smithsonian Report of 
1871 ; pp. 326 and 335.) It can be shown that the adjustment formulas 
thus constructed are embraced under the general form 
[(2m + 1)K) — 2(m 2 )(m 2 )>' = (m 4 )w +[(ro 4 ) — l 2 (m 2 )]K + M -i) 

+[(»*) — 2V)](%+«- 2 )+ • • • • 

+[(m 4 ) — m%m 2 )Ju m +u_ m ) (3) 

where the notation used is 

(m») = 1" + 2" + + m". 
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The values of the coefficients I for all values of 2m + 1 from 5 to 25 inclu- 
sive, are here given in Table A. They are only carried out to two or three 
places of decimals, which are all that are needed in practice, and the last 
decimal figure is in some cases altered by a single unit in order that the 
condition / + 2^ + h -+■■••• + 40 — 1 may be exactly satisfied. 

Table A. 

(giving minimum values to e'~-e.) 



2m+l 


5 


7 


9 


11 


13 


15 


17 


19 


21 


23 


25 


h 

h 

h 

h 

h 

h 

h 

h 

h 

ho 

hi 

'l2 


.48 

.34 

-.08 


.34 

.29 

.14 

-.10 


.26 
.23 
.17 
.06 
-.09 


.208 
.196 
.161 
.102 
.021 
-.084 


.174 
.168 
.147 
.112 
.063 
.000 
-.077 


.152 
.147 
.133 
.110 
.079 
.038 
-.012 
-.071 


.134 
.130 
.121 
.105 
.083 
.056 
.022 
—.019 
—.065 


.118 
.117 
.110 
.099 
.084 
.064 
.039 
.011 
—.023 
—.060 


.108 
.106 
.101 
.093 
.081 
.067 
.049 
.027 
.003 
—.025 
—.056 


.098 

.097 

.093 

.087 

.078 

.067 

.053 

.037 

.019 

—.002 

—.026 

—.052 


.090 

.089 

.086 

.082 

.075 

.066 

.056 

.043 

.028 

.012 

—.006 

—.027 

—.049 


e 


.69 

.29 


.58 
.19 


.50 
■13 


.456 
.109 


.418 
.0907 


.389 
.0782 


.365 


.345 


.328 


.313 


.300 


K 4 ) 


.0677 


.0600 


.0548 


.0496 


.0456 



The probable error of a term adjusted by any such formula as (1), assu- 
ming that the true law of the series is algebraic and of a degree not higher 
than the third, will be 

e' = i/PS«3+*i(«!+£i)+ • • • +&(4+«i,)] 
where e , e 1 &c, denote the probable errors of the given terms u , u 1 &c. 
If these are all equal, denoting their common value by e we shall have 



£' 

e 



^jQg + 2® + % 



+ 



+ l\ 



o], 



(4) 



and this may be regarded as approximately the ratio which the probable 
error of any adjusted term will bear to that of the corresponding unadjusted 
one, so far as this ratio can be determined a priori. It is true that the er- 
rors of the given terms will not in general be really equal, but it usually 
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happens in statistical and other series that the weights of the terms follow 
a tolerably regular sequence, so that if for instance the errors £j, s 2 &c., are 
greater than s , the errors s_ x and e_ 3 will be less than e , and vice versa, 
and thus they tend to correct each other. The values of I 2 too are much 
larger for the middle of a formula than for its extreme terms, so that the 
error s' is chiefly due to the errors of u and a few other terms adjacent to 
it. The values of e'-e- e corresponding to the several values of 2m + 1 in 
Table A. are given at the foot of that table. These are minimum values, 
for the reason that, under the conditions we have assumed, the method of 
least squares gives the most probable value of the adjusted term u' . 

Before having met with Schiaparelli's work the present writer had 
constructed, with special reference to the adjustment of mortality tables, 
various systems of formulas of this general character, one of the most ad- 
vantageous of which seemed to be that which renders the probable values 
of the fourth differences of the adjusted series a minimum. (Smithsonian 
Report of 1871, p. 332.) The series is thus graduated with the greatest 
possible smoothness. The coefficients under this system are here given in 
Table B. In a continuation of the article just cited, (Sm. Rep. of 1873, p. 
327), the subject was still farther investigated. The fourth difference of any 
five consecutive terms in the given series being 

Table B. 

[giving minimum values to (J' 4 )-4-(J 4 .)] 



2m+l 


5 


7 


9 


11 


13 


15 


17 


19 


21 


23 


25 


h 

h 

h 

h 

h 

h 

h 

h 

h 

h 

no 

hi 

'l2 


.56 

.29 
-.07 


.42 

.29 

.05 

-.05 


.34 

.27 

.11 

-.02 

-.03 


.290 
.245 
.135 
.026 
-.028 
-.023 


.252 

.222 

.145 

.056 

-.007 

-.027 

-.015 


.222 

.202 

.146 

.076 

.015 

-.018 

-.022 

-.010 


.200 

.184 

.143 

.087 

.034 

—.004 

-.020 

-.017 

—.007 


.182 

.170 

.137 

.093 

.047 

.010 

—.012 

—.018 

—.013 

—.005 


.166 

.157 

.132 

.096 

.057 

.022 

—.003 

—.015 

—.016 

—.010 

—.003 


.1534 

.1463 

.1262 

.0969 

.0635 

.0317 

.0065 

-.0092 

-.0151 

-.0133 

-.0077 

-.0025 


.1424 

.1368 

.1205 

.0962 

.0678 

.0395 

.0151 

-.0024 

-.0119 

-.0140 

-.0110 

-.0060 

-.0018 


e 


.70 
.23 


.60 
.080 


.54 


.495 


.461 
.0104 


.433 
.0064 


.409 


.389 
.00273 


.372 
.00185 


.358 


.345 
.00091 


(A) 


.038 


.0187 


.00422 


.00124 
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J 4 = 6w — 4(« 1 +«_ 1 ) + (w 2 +m_ 2 ) 
its probable value, without regard to sign, is 

(J 4 ) = e ,/[62 + 2(4 2 +l)] = Sl /70, 

and hence we have 

e = .11952(J 4 ) (5) 

It was shown (p. 330) that the probable value of the fourth differences 
of the adjusted series is 

(J' 4 ) = Sl /M+2(4i+4H....+4«)] ) • • • (6) 
where J Q , J 1} d 2 &c, are the fourth differences of the series 

^2> 'l' 'o> I'D '"a t 3 . . . . t m , 0, 0, 0, U. 
Combining (5) and (6), we get 

(^4 = .11952 1 /[^ +2(J?+J|+....+^ +2 )], . . .(7) 

which enables us to compute, for any adjustment formula such as (1), the 
ratio which the probable value of any fourth difference of the adjusted se- 
ries bears to that of the corresponding fourth difference in the unadjusted 
series. The values of this ratio for Tables A. and B. are found at the 
bottom of the tables. They measure the smoothness of graduation. For 
instance, the ratio .00091 for the 25-term formula of Table B. shows that 
when that formula is used the adjusted series will be smoothly graduated 
to three places of decimals farther than the given series is. 

An error was committed in connection with this subject, which had es- 
caped my notice until quite recently, and I take this earliest opportunity to 
correct it. Since we have for the given series the relation 

s = .11952(i 4 ) 
it was too hastily inferred that we shall also have for the adjusted series 

s'=.11952(J' 4 ) 
and consequently 

g' _(A) 
« (^) 

from which the conclusion was drawn that the coefficients I in Table B. 
(called Table III. in the article referred to) must render the ratio e' -*- e a 
minimum, in as much as they are known to render (J' 4 )-h( J 4 ) a minimum. 
It was thus made to appear that the formulas of Table B. are the ones 
which give the most probable value to the adjusted term u' , for any given 
value of 2m -f 1, and that the ratio (z/' 4 )h~(z/ 4 ) is, in theory at least, equal 
to the ratio of the probable error of the adjusted term to that of the given 
one. A great discrepancy between this result and the facts of the case was 
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indeed recognized at the time and remarked upon in the context, but it was 
accounted for by the existence of causes which were really of minor impor- 
tance, the true cause, an error of theory, being overlooked. The truth is 
that, according to the theory of probabilities, the process by which we ob- 
tained the relation 

e = .11952(J 4 ) 
implies that the errors of the terms u , u l , &c, of the given series are 
wholly independent of each other, and that errors in excess and errors in 
defect are equally likely to occur. This is really true for the terms of the 
given series, but it is not true for the adjusted series, because any five con- 
secutive adjusted terms have been all computed, to a considerable extent, 
from the same given terms, only taken in somewhat different proportions, 
so that the errors of these adjusted terms will most probably be all in the 
same direction, that is, all in excess or all in defect. Hence the relation 

af = .11952(J' 4 ) 
does not hold good, even approximately, and any deductions from it must 
fall to the ground. The formulas given in pages 330 and 331 of the Smith- 
sonian Report of 1873 will be corrected if we substitute (J' 4 )-r-(z/ 4 ) for s'-s-s 
where ever the latter occurs. So also the value 0.232 on page 332 really 
represents the ratio (J' 6 )-h(J 6 ), and the values 0.305 and 0.460 on page 333 
represent the ratio (^ / 2+2)~H^2+2)- They measure the smoothness of the 
adjustment, but they do not measure its accuracy. And it is not true as 
stated at the foot of p. 337, that when repeated adjustments are made^ the 
ratio of the probable error of the final series to that of the original one 
will be the product of the ratios due to the two or more formulas used. To 
find the probable error in such a case, we shall have to express the final ad- 
justed term as a linear function of the original terms, and then proceed as 
in the case of formulas (1) and (4). 

The correct values of the ratio of error e'-r-e, as found by formula (4) 
from the coefficients in Table B., are given at the foot of that Table. They 
are somewhat larger than those of Table A., while the ratios of irregularity 
(^' 4 )-i-(J 4 ) are a good deal smaller. It is impossible to get at once the 
maximum of probable accuracy and the maximum of smoothness. Formu- 
las may, however, be constructed so as to give intermediate and moderately 
small values to both ratios. Very good results are obtained by assuming 
that the coefficients I are ordinates to the curve 

the constant interval between these ordinates being taken as the unit of 
x. The values of the four constants are found from the four equations 
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(8) 



(2m+l)^+2[(m<)C+(m 6 )Z>+(m 8 ).E] = 1 "^ 
(m 2 )J.+(m 6 )C , +(m 8 )DH-(m 10 )£ =0 

2C+ 3(m+l) 2 D+4(m + 1) 4 ^ = _ 

The first two equations follow from Schiaparelli's conditions (2) already 
given, while the last two require that y and dy~-dx should both become zero 
when a; = ± (m + 1). (Compare the 8m. Beports of 1871 and 1873, pp. 
322 and 352.) Having obtained the numerical values of A, C, D and E 
for any assumed value of 2m -f- 1, and substituted them in the equation of 
the curve, we assign to x the values 0, 1, 2, ... . to, in succession, and the 
resulting values of y will be l , l x , l 2> . . . . l m . These are here givea in 
Table C, for values of 2m+l from 5 to 25. The ratios s'-t-e are but little 
larger than those in Table A., while the ratios (J' 4 )-f-(J 4 ) are considerably 
smaller, excepting for the very shortest formulas, which are of no great con- 
sequence. 

Table C. 



[giving intermediate values to s'-f-e and {4' i )-r-(J i ).'] 



2m-f-l 


5 


7 


9 


11 


13 


15 


17 


19 


21 


23 


25 


*n 


.44 


.32 


.24 


.206 


.176 


.154 


.136 


.122 


.112 


.102 


.094 


'i 


.37 


.30 


.24 


.204 


.176 


.154 


.136 


.122 


.112 


.102 


.094 


U 


-.09 


.14 


.18 


.177 


.162 


.146 


.133 


.121 


.110 


.102 


.094 


h 




-.10 


.03 


.093 


.117 


.121 


.118 


.111 


.104 


.098 


.091 


h 






-.07 


-.020 


.038 


.070 


.085 


.090 


.090 


.088 


.085 


h 








-.057 


-.037 


.005 


.037 


.056 


.067 


.071 


.073 


h 










-.044 


-.040 


—.012 


.015 


.034 


.047 


.054 


I, 












-.033 


—.039 


—.021 


.000 


.018 


.031 


h 














—.026 


—.035 


—.025 


—.009 


.007 


h 
















—.020 


—.032 


—.027 


—.015 


no 


















—.016 


—.028 


--.026 


'n 




















—.013 


—.024 


l\i 






















—.011 


e' 

e 


.70 
.36 


.58 
.19 


.50 
.088 


.461 
.0516 


.427 
.0324 


.398 
.0207 


.375 


.354 


.339 


.325 


.311 


(4'*) 


.0146 


.0103 


.00802 


.00589 


.00480 


{4 a) 

























The question as to which of the formulas of the three tables A, B or C 
is most advantageous for use in any given case, would be easily answered if 
the conditions we have assumed were really fulfilled in practice, that is to 
say, if the series of statistical or physical observations which we wish to 
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adjust were really of algebraic form and of a degree not higher than the 
third, subject only to accidental deviations from this law. The best formula 
would then be the 25-term formula of Table A., since this has the smallest 
ratio of error e'-f-e. But actually, the true law of an observed series will 
probably be quite unknown, and we can only assume that a group of 2m+l 
terms will follow the algebraic law approximately, and this assumption may 
not be safe for so many as 25 consecutive terms. The effect of a change in 
the form of the function, or law of the series, may be illustrated thus : 
Take the following equidistant values of the function 



u = sin i 



which is not of algebraic form. 



e 


u 


e 


M 


e 


u 


15° 


.25882 


35° 


.57358 


60° 


.86603 I 


20 


.34202 


40 


.64279 


65 


.90631 


25 


.42262 


45 


.70711 


70 


.93969 


30 


.50000 


50 
55 


.76604 
.81915 


75 


.96593 



Apply to these values the 11-term formula of Table A. and the 13-term 
formula of Table B. The two adjusted values of the middle term sin 45° 
will be .70702 and .70712 which differ from the true value by .00009 and 
.00001. 

The formula from Table B. gives the best result, although its error-ratio 
s'-r-e is a little the largest. The reason is mainly that l and l t are larger, 
and l m is smaller, in that formula than in the other, giving a greater weight 
to the middle term and terms adjacent to it. Thus the formulas of Table 
B. follow the actual curvature of the given series more closely than those of 
Table A. It seems advisable not to be governed solely by the theoretical ratios 
s'—e, but to employ such formulas as appear advantageous on various grounds, 
and rest content if the series thus adjustsd is a smoothly graduated one, and 
satisfies certain tests of good adjustment, that is, certain conditions which 
would probably be satisfied by the true series, if we had it. There is some 
advantage in having a series graduated smoothly, even if the probable er- 
rors are such that two or three of the last decimal figures are in doubt, and 
cannot be regarded as having strictly any real value. The chief advantage 
is, that errors of computation can be more readily detected by the irregu- 
larity they produce in the differences of the series. 

(To be concluded in No, 4.) 



